Spot deformation and replication in the two-dimensional Belousov-Zhabotinski 

reaction in water-in-oil microemulsion 
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In the limit of large diffusivity ratio, spot-like solutions in the two-dimensional Belousov- 
Zhabotinski reaction in water-in-oil microemulsion are studied. It is shown analytically that such 
spots undergo an instability as the diffusivity ratio is decreased. An instability threshold is derived. 
For spots of small radius, it is shown that this instability leads to a spot splitting into precisely 
two spots. For larger spots, it leads to deformation, fingering patterns and space-filling curves. 
Numerical simulations are shown to be in close agreement with the analytical predictions. 
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Localized patterns such as spots belong to the class 
of dissipative structures found far from equilibrium 
In recent years, considerable progress has been made in 
the understanding of these systems. The question of sta- 
bility of such patterns is central and the source of in- 
stabilities must be carefully scrutinized. In particular, 
the occurrence of instability can lead to the deformation 
of spots followed by spot multiplication (also called self- 
replication) or fingering. This intriguing phenomenon 
has been the subject of research since the pioneering work 
of Pearson Q. Shortly after, thanks to the development 
of open spatial chemical reactors, self-replication was 
observed in various experiments such as ferrocyanide- 
iodate-sulphite reaction the Belousov-Zhabotinsky 
reaction [J, |5|, |6| , and chloride dioxide-malonic-acid re- 
action Q. By now this phenomenon is believed to be 
universal [1, @| . It is not restricted to chemical reactions 
and occurs in many systems in bio logy [Io| , material sci- 
ence [III . [HI and nonlinear optics |l3l |. 

Analytically, spot replication is relatively well under- 
stood in one dimensional setting. Nishiura and Ueyama 
[l4j ] proposed that self-replication of spikes occurs when 
the spike solution dissapears due to the presence of a 
fold point. Similar explanation has been reported for 
the box-like patterns [15|. In two dimensions, a mech- 
anism for spot instability has been proposed in Q for 
general reaction-diffusion systems; related analysis was 
performed earlier in the framework of a piecewise-linear 
approximation in [l6| , [ItJ and more recently in the con- 
text of diblock copolymer systems 11 1 , [12] • See also [l|| , 
[2fj| and 21 1 and for a related approach from the point 
of view of interface motion. 

In this letter we perform an analytical and numeri- 
cal investigation of the two-dimensional localized spots 
that were recently reported for the Belousov-Zhabotinski 
(BZ) reaction in water-in-oil microemulsion [Bj], [||. In 
the classical BZ reaction, spiral waves are observed [HI 
but no localized spot solutions are possible. Indeed, lo- 
calized structures develop only when the ratio of diffusion 



coefficients is sufficiently large, which occurs in the mi- 
croemulsion system but not in the classical BZ reaction. 

We consider the water-in-oil microemulsion model of 
the BZ reaction as described in Q and 
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where v, z are dimensionless concentrations of activa- 
tor HBrC>2 and oxidized catalyst [Ru(bpy)3) 3+ respec- 
tively; D v and D z are dimensionless diffusion coeffi- 
cients of activator and catalyst; /, e and q are parameters 
of the standard Keener-Tyson model [l8[; io represents 
the photoinduced production of inhibitor, and m repre- 
sents the strength of oxidized state of the catalyst with 
< mz < 1. This reaction was shown experimentally 
and numerically to admit localized spot patterns that 
persist for long time [B| , @ . 
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We rescale the variables as z — 1 jm - 
to -1 / 2 ?), t = e^m 1 ! 2 !. In the new variables, after drop- 
ping the hats, we obtain 

v t = e 2 Av + f(v, w); rw t — DAw + g(v,w) (2a) 

where 



f(v, z) = - [f + fiw 
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and with the nondimcnsional constants given by 
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s D v m 1/2 , D = D z e\m 
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In the limit m — > oo, /i — > and r — > we obtain the 
reduced system, 

vt = e 2 Av — fo + wv — v 2 

v + q (4) 

= DAw + l-vw. 

In particular, parameter values used in Fig. 14 of @ 
are fo = 2.18, io = 0, m = 10, ei = 0.01, eo = 0.1 and 
D z /D v = 100 which gives a = 0.3, /i = -0.007, r = 
0.01, so that the simplification (jU is appropriate. 

Experimental and numerical evidence in || and Q 
suggest that ([T|) admits localized spot soltuions, such 
as shown in Fig. [T] Such spots occur in the regime 
where £ « 1. Our goal is to describe analytically the 
radius and profile of such a spot and then study its sta- 
bility using singular perturbation techniques similar to 
those described in Q. As we will demonstrate, the in- 
stability thresholds appear in the regime where D ^> 1. 
For sufficiently large values of D a spot pattern is sta- 
ble. However as D is decreased, instabilities of the form 
exp(Ai) cos(m8)(f)(r) may develop, where 9 and r are the 
angular and radial coordinates, respectively. We provide 
the analytic description of the profile of the spot (0 fTOj) 
and the dispersion relation (Ti"4|) between m and A. This 
relation leads directly to the estimate for the instability 
threshold. 

We begin by constructing a stationary (time- 
independent) spot-type solution on a two dimensional 
unit disk {x : \x\ < 1}. We assume that D ^> 1. Then 
to leading order, w ~ wq is a constant to be determined, 
and the solution for v consists of an interface located at 
| a; | ~ I which connects two nearly spatially homogeneous 
layers. To find the profile of such an interface and its 
location I, let us rescale near I as v(x) = V (y) with 
y — (r — I) /e and r = \x\ . For the steady state we then 
obtain, to leading order, V" (y) + f (V; w ) — 0. The 
interface solution corresponds to a heteroclinic orbit of 
this ODE. The existence of such an orbit is only possible 
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FIG. 1: (a) Radially symmetric, stable spot solution of the 
reduced BZ system Q in the BZ system on a unit disk, (b) 
Radial profile of v (solid curve) and its asymptotic approxima- 
tion (dashed curves) given by (0. Insert: a blowup showing 
the profile of the interface, (c) Radial profile of w (solid curve) 
and its two-term asympotic approximations (dashed curves) 
given by (|10|l . Parameter values are D = 20, fo = 2.0, e = 
0.02, q = 0.005. 



whenever the equations 
f(V-:z )=Q = f(V+:w ) 



f{ S] wo)ds = (5) 



are simultaneously satisfied for some values V- < V+. 
Since q <C 1, we have V- ~ and (O can be written as 



(6) 



Next, we ignore the 0(q) terms and integrate V"(y) 
f(V,Wo) =0 to obtain 
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with Kf. and wo given by (|t>J) ana wliere r — \x\ 
formula describes the profile of the interface in the limit 

— 1/2 

q — * 0. Its thickness is of 0(eV + ). To determine its 
location I, we integrate the second equation in (|2a|) . Zero- 
flux conditions then yield / g = so that 

g (V + ,tv )l 2 + (l~l 2 ) g(0,w o )~0 (8) 

In the limit of the reduced system ([4]), we obtain 
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To determine the correction to w, we write w 
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D~ 1 wi + O (£> ) • to obtain Awi + g (w , v ) = 0. Im- 
posing continuity at the interface r = I, the solvability 
condition wi(l) = 0, and using ([5]) and g (0, Wq) = 1, we 
then obtain 

d-W-r!) r<l 

■ (10) 
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An example of a localized spot and its radial profile are 
shown in Fig[TJ 

When decreasing the diffusion coefficient D of the re- 
covery variable, numerical simulations show that the spot 
becomes unstable. To compute the threshold associated 
with this instability, we linearlize around the localized 
spot solution ([7]) as 

v{x, i) — v(x) + exp (At) cos (mff) 4> (r) 
w(x, t) — w(x) + exp (At) cos (m9) i\) (r) 

where v and w are given to leading order in ([7] and [9]) 
spot-type solution, m is an integer, ^ < 1 and (r, 9) 
are the polar coordinates of x. Substituting into |T]) we 
then obtain 
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Note that v r satsifies that e 2 (v rrr + l/rv rr — l/r 2 v r ) + 
f v v r + w r f w = 0. Therefore multiplying (jllbp by v r and 
integrating by parts, we obtain 



A 



(m 2 — l)e 2 



Vrfw (tp ~ W r ) ■ (12) 



where we have assumed that to leading order, 
4> ~ v r ,\ <C l,ip <C 4>- Since w r is ex- 
ponentially small outside the interface, we simplify 

/ V r fw (V 1 ~Wr) ~ - (0 - W r (I)) f w (v, W ) dv . 

We estimate w r (l) ~ — l/(2D)g (V+, wq) and to deter- 
mine ip (I) , we integrate (jllbp over the interface I to ob- 
tain Dip r \\- ~ — <? (0, wo) /i 2 where we have used ~ v r 
and ([H]). Keeping only leading order terms in D we then 
obtain the following problem for 



= ?/Vr + -tpr o"^ = 0, 



(13a) 



^(o) = o = #(i), v' (^ + ) - V-' (r) = (o, wo) 



(13b) 



Using the continuity of ip at r — I we get 
1 



ff (0,w )(l + P™) 
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and substituting into (j!2p we obtain 
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2eD efv? 
For the reduced model (j4j) we obtain an explicit result 
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In Fig. [21 the analytical prediction given by (|14p is found 
to be in good agreement with the numerical computations 
of problem (fTTj) . Full numerical simulations of Q also 
agree with this prediction. For example, when taking 
e = 0.03, fo = 1.3 and q = 0.01, slight spot deformation 
corresponding to mode m = 2 is observed when D = 1.0 
but not when D — 1.3. This agrees well with D ~ 1.1, 
the threshold predicted by (fT4")l . 

From (fl~4|) and (|16p . it is clear that the mode m = 1 
is always stable and that for large enough £>, all modes 
m > 1 are stable. As 13 is decreased, instability sets 
in when D = 0(e^ 1 ). The threshold value is found by 
simultaneously solving A = dX/dm = 0. The resulting 
graph is shown on Fig. [^Ji. In particular, note that the 
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FIG. 2: (a), (b) Comparison of numerical computations of 
A given by l|lip (diamonds) with the analytical result (1141 
I16[) (dashed line) for the reduced model j|} on a unit disk. 
Parameter values are q = 0.01, e = 0.01 and (a) D = 3, fo = 
1.3; (b) m = 2, fo = 1.3. (c) Simultaneous solution of A = 
= dX/dm, showing the first value of D for which instability 
occurs. The system is stable above the curve and unstable 
below it. To compute A numerically, (|ll|l was reformulated as 
a boundary value problem by adjoining the equation dX/dr = 
along with fixing tp(l). Maple's numerical boundary value 
problem solver was then used with initial guesses </> — v r , X = 
and ip = the solution of (|13p . All computations are correct 
to four significant digits. 



system is stable if De > 0.038, independent of the value 
of fo- In the limit of small radius I — > 0, the first unsta- 
ble integer mode is to = 2 corresponding to A = 6, so 
that the spot of small radius becomes unstable whenever 
eD/ 3 / 4 < 0.0727. More generally, by eliminating A we 
find that there exist constants l± < 1% < . . . < 1 such the 
first unstable mode is to provided that l m -i < I < Im, 
where h = 0, l 2 = 0.491, l 3 = 0.667, h = 0.753 and 
lm ~ 1 — 0.937/to as to — > oo; where 0.937 is the root of 

(14) e- 22 (3 + 2z) + 3-4z = 0. 

Numerical computations indicate that self-replication 
is more prevalent for spots of small radius (see Fig. [3]). 
For larger spots, a deformation usually leads to to the 
so-called "finger growth" and space filling curves. Others 

(15) studies have shown the occurrence of fingering instabil- 
ities leading to labyrinthine patterns [22] , [HI, The 
question of whether the self-replication or fingering in- 
stability occurs first is still open. 

For smaller, O(l) values of D, there is also a differ- 
ent instability mechanism that can lead to splitting of a 
spot into a ring as illustrated in Fig. [3J Unlike spot mul- 
tiplication, this instability is radially symmetric, and is 
caused by the dissapearence of the steady state solution 
- whereby the steady state solution ceases to exist due to 
the presence of a saddle-node bifurcation - rather than 
by its lateral instability. Numerical simulations suggest 
that that spot replication occurs only for spots of smaller 
radius, whereas spot-to-ring instability is dominant for 
larger spots. 

To conclude, we have estimated analytically for which 
value of diffusion D spot deformation first occurs. As 
is D is decreased further, self-replication and/or spot-to- 
ring instability is observed. As evidenced by numerical 
simulations, we conjecture that spot deformation is the 
precursor to this phenomenon. 
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FIG. 3: (a) Bifurcation diagram for Q in D and fo with e = 0.05, q = 0.01. Solid dots represent deformations of a spot without 
topological change. Points marked by "2" represent spot-replication into two spots. An empty circle represents spot-to-ring 
bifurcation and empty circle with a number inside represents spot-to-ring-to-spots bifurcation. A solid line represents the 
boundary of spot-to-ring replication which occurs at the fold point of the radially symmetric steady state. The bifurcation 
diagram was obtained by solving the full two-dimensional system Q using the finite element package FlexPDE [23| with zero- 
flux boundary conditions and 800 elements on a quarter-disk. For initial conditions, {7J |9] 1 10 p was used, but with r replaced 
by r(l + 0.05 cos 28). The solid line was obtained by solving for the fold point of the radially symmetric steady state using 
Maple's boundary value problem solver, (b) Full numerical simulations starting from a spot-state (0 1 10 p as initial conditions. 
First row: spot deformation, fo = 1.3, D = 0.35. Note that the final stedy state of the system is a deformed blob shown here 
at t = 4000. Second row: self-replication, fo — 2.2, D — 0.24. Third row: spot-to-ring bifurcation, fo = 1.3, D = 0.1. Fourth 
row: spot-to-ring-to-spots bifurcation, fo = 2.6, D — 0.1. 
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